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Abstract 

The spectrum of collective fermionic excitations in a finite temperature 



J> ■ QED34.1 is studied in different regimes. It is shown that within the standard 

00 _ 

^ ■ perturbation approach the one-loop dispersion equation, besides the ordinary 

one-particle excitation, has four new solutions. The additional excitations are 

00 

ON ■ gauge-dependent and two of them have nonphysical signs of residues in the 

propagator poles. The temperature evolution of the solutions is investigated 
and it is shown that the use of effective propagators leaves no more than one 
additional mode which becomes propagating at WM, when the gauge 
invariance is restored. The other three modes, including those with nonphys- 
ical residues in the propagator poles, are always strongly damped, thus the 
thermal effects do not produce pathologies in QED3 + i. 
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I. INTRODUCTION 



A significant advacement in studying gauge theories at high temperatures has been 
achieved in recent years (see e.g. Ref. []]]). Hard thermal loops (HTL) resummations tech- 
nique has made possible to penetrate into the temperature region, where the effective 
coupling constant is large, gT/M 3> 1 (M is the vacuum fermion mass), and the ordinary 
perturbative series fail. The low-temperature limit, where the expansion in the vacuum 
coupling constant is valid, is also well-known (see, e.g. Ref. [[|). However, the intermediate 
range of temperatures comes as surprise: the naive one-loop calculation in the framework of 
QED shows that at T ^ 4M new modes arise in the fermion spectrum. 

At low temperature the modification of the fermion dispersion equation is obvious: ther- 
mal effects produce only a small shift of the vacuum mass shell of the one-particle excitation 
(OPE) and no new branches in the spectrum are formed. At high temperature the loop cor- 
rections may be of a tree-level order. For instance, when calculating the fermion self-energy 
every perturbation order brings the factor T 2 (which is typical for diagrams containing hard 
thermal loops, electron self-energy among them) and the effective coupling constant gT/M 
becomes (unlimitedly) large and one must sum out HTL. In practice, HTL summation im- 
plies substitution of the effective propagators, which are functions of the effective masses 
and dampings instead of the bare ones. In (extremely) high temperature limits, when the 
vacuum fermion mass is negligible, an additional hole-like state with a negative frequency 
arises in the electron branch, as well as a state with a positive frequency in the positron 
branch. Thus the number of excitations doubles as compared to the low-temperature case. 
The additional solutions (collective excitations, CE) were found in Refs. where the 

one-loop self-energy of (massless) fermion was calculated (the spectrum structure at high 
temperatures has been studied also in Refs. [0|§). Masses of the above excitations are of the 
order of gT, and in the leading approximations the corresponding solutions to the dispersion 
equation are gauge-invariant 



Early in the 90th it was shown in Ref. M in the high-temperature limit and in Ref. [1C 
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at the "naive" one- loop level in QED with massive fermions that at T > AM in each 
branch a new low-frequency mode arises. Besides, it was demonstrated in Ref. |J that the 
above excitation disappears after the damping of OPE is taken into account in the effective 
propagator. 

Therefore, several questions arise: how high must the temperature be to let collective 
excitations come into existence, how do they exactly evolve with temperature, and are 
they physical excitations or not? In the present paper we shall study properties of the 
additional solutions to the dispersion equation, point out the physical ones, and estimate 
the temperature required to make them propagating. Within the naive perturbation theory 
the additional modes arise at rather low temperatures, T «1 AM. These modes are situated 
in the low frequency range (oj ~ M 3 /T 2 at zero momentum) and always arise in couples: two 
modes at positive and two at negative frequencies, both on electron and positron branches. 
However, the additional solutions are pathological: one excitation in each pair has the 
"wrong" sign of the residue. Besides, all the modes are (in leading order) gauge dependent. 
When the temperature is raised the additional positive-frequency modes (in the electron 
branch) and one of the negative-frequency modes (that with the wrong residue sign) remain 
in the low-frequency region, while the absolute value of frequency of the other modes rise 
and becomes ~ gT. 

Since the difficulties described above jeopardize not an abstract model, but quantum 
electrodynamics, the theory describing the real world, one should believe the pathological 
modes (those which are gauge- variant and/or have indefinite metrics) should be exorcised 
by applying a correct calculation procedure. 

We shall demonstrate that after inserting the damping of OPE into the bare propagators 
among positive frequency exitatons only one mode survives at all temperatures. At T ^ 5M 
the real part of the inverse propagator has two zeroes at negative frequencies, but its imagi- 
nary part is larger than the frequency, thus the physical excitations are absent, and only at 
10M one of the negative-frequency solutions overcomes the damping and becomes the 
propagating one. This is in a good agreement with the approach to the region T ~ 10M 
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from the higher temperature domain where the hard thermal loops results are valid. 

The paper is organized as follows. In Sec. 2 the one-loop self-energy in an arbitrary 
relativistic gauge is calculated using the bare propagators. We show that the above naive 
calculation brings pathological modes into the model. In Sec. 3 we study how the damping 
may affect the collective excitations in QED plasma and show that the "improved" naive 
approach and HTL scheme are in a good agreement with one another in the intermediate 
temperature range. A possible physical interpretation of the additional high-temperature 
excitation found in the fermion spectrum is given in Sec. 4. 



II. ONE LOOP ELECTRON SELF ENERGY AT FINITE TEMPERATURE 

In this section we shall study the spectrum of one-particle fermion excitations within the 
standard one-loop approximation, i.e. using the free electron and photon propagators. We 
shall use the thermal Green functions method PJTT|-|i"3"! (of course, the real-time approach 
14] , [T5| JT|| produces the same results). The spectrum is defined as poles of the Fourier- 



transformed retarded Green function Sr(p , p) = Sr(p) or, equivalently, as zeroes of S R ^ 1 (p). 
The propagator Sr(po,p) is evaluated from the thermal Green function S(iu> n ,p), uj u = 
(2n + l)7r/3, (3 = 1/T via analytical continuation iu n — > uj = po, holomorphic in the upper 
half-plane of uj and having correct asymptotic behavior at uj — > 00 (the fermion self-energy 
falls as 1/uj at high frequencies) PJlB] . 



The Lagrangian of QED 3+1 in a relativistic gauge a is 

C = --F^ --^-(d^f + ^S + gA-M)^ . (1) 
4 la 

The fermion self-energy E(p) is defined as []: 

s-\ p ) = s 1 (p) + j:(p), po = iu n , (2) 



1 We define the Green functions (propagators) in the temperature technique as average of the 
T-products in thermal time r of field operators. 
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In the one-loop approximation it equal to 



r rlk 9irl 

£i(p) = -! S 2 T £ /(|^V?„(j> + k) lu D"(k) , i„ = i^ , (4) 

The study of £i(p) for arbitrary temperature, frequency and external 3-momentum p 
is quite difficult and we consider the case of zero spatial momentum, p = 0, only The 
corresponding expression for Si may be written as (in Ref. [JTOj] £i(co>) was calculated in the 
Feynman gauge): 



S 1 (^,0)^S 1 H=7°(aH + (l-a)6H) + ^cM , (6) 
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V-|„ _ x \ • CD 



g 2 u dk (k{M 2 - u 2 + 2k 2 ) 2k 2 e \ 

a ~^ 2 ~Jo Au 2 k 2 - (M 2 - uj 2 ) 2 { e? k - 1 + eP e + 1 J ' 



g 2 uj(M 2 -u 2 ) y°° * 

2tt 2 7o (4cj 2 A; 2 - (M 2 - c^ 2 ) 2 ) 2 lJ 



'£;(8^ 2 £; 2 - 4A; 2 M 2 - (M 2 - c^ 2 ) 2 ) P((M 2 + c^ 2 ) 2 - 4e 2 M 2 ) 



\ e /3fc _ ! e f3e + ! y ' 

_2g 2 M dk, ( k{M 2 - u 2 ) k 2 {M 2 + u 2 )\ 

°~~^ 2 ~lo Auo 2 k 2 - (M 2 - uo 2 ) 2 \ e? k - 1 + e(e^ e + 1) J ' l0) 

the vacuum contribution, which is insufficient at high temperatures, is omitted). 
The spectrum is defined as nontrivial solutions to the equation 

( 7 °o;-M-£i(u;))V> = . (11) 



We suppose that an excitation belongs to the electron branch if if) may be written as 



u 



1>=\ Q \ , (12) 



and to the positron branch if if) is 



H°) • (i3) 

We shall consider below the electron branch excitations only, since the positron branch 
solutions may be obtained by reversing the sign of u. The electron spectrum is determined 
as solutions to the following equation: 

u — M = a(u) + (1 - a)b(u) + ^c(w) . (14) 

In the vicinity of a pole the retarded propagator may be presented as 

S R (w) ~ Z E — — — , 15 

uo — E + io 



to = E is a solution to the dispersion equation, ipE is the solution to the Eq.flTTD for uo = E 
normalized to unity, and the residue Z E is 

1 d 



(16) 

uj=E 



At low temperature corrections to the free propagator are small and the dispersion equa- 
tion has just one solution, uo w M. When the temperature is raised, Hx(oj) increases as T 2 , 
but for low frequencies E%(uj) varies drastically due to the contribution to the integral of the 
momentum region corresponding to pole in the integrand (this happens since the poles in 
the photon and electron propagators in the expression for S 1; Eq.([|) are very close to one 
another). 

In Fig. 1 the dependence of real part of the fermion self-energy on frequency for the 
electron branch in the Feynman gauge a = 1 (i.e. the function a(u) + c(u))) is presented for 
T = 3M and 5M. The intersections of 3?eSi with the straight line uo — M are solutions to 
the dispersion equation in the electron sector. One can see that starting with T ^ 4M there 



are five solutions in each branch. It follows from the continuity of S(u;) as a function of to 
that the additional solutions come by two (the plot of ^teLi at T = 5M in Fig.l confirms 
this) and that the residues of two poles in each pair have opposite signs. Actually, the exact 
location and the residue of the poles depend on the gauge. This follows from the explicit 
form of the function Si (a;) Eq.(|B|). Consider the pole closest to zero as an illustration. Its 
location and residue may be obtained by expanding Si(cj) in series in u (the coefficients are 
calculated in the leading order in T/M): 

M + 

w = T7(to • * = i + c-«>9 • m 



with 



' k 3 dk (—- + ——) « -a'(0) « -6'(0) , (18) 



vr 2 M 4 Jo Ve fc - 1 e fc + 1, 

2^2 T 2 f°° , , f 1 1 



«i, ** to + —i) " c <°> • < 19 > 

i.e. parameters of this excitation do depend on a. 

Thus we see that the additional excitations arising in the perturbation theory are patho- 
logical: their characteristics are gauge-variant and some residues are negative, which implies 
indefinite metrics. One should believe that when a correct calculational procedure is used 
the excitations with unacceptable properties would vanish from the physical spectrum. 

When the temperature increases excitations behave in a differentl way. Energies of the 
middle three remain small. This may be easily understood: the above excitations arise 
due to large oscillations of Si in the small frequency region. As we have discussed above, 
those oscillations are due to large contribution of the momenta corresponding to the pole 
of the integrand in Si. It is well-known that at high temperatures the effective mass of 
OPE becomes much greater than the vacuum one. If one neglects the vacuum mass while 
calculating S 1; the dispersion equation has two solutions for the electron branch (and two 
for the positron one), then three middle excitations vanish, and the energies of the remaining 
modes grow when the temperature is raised and the corresponding poles prove to be gauge- 
invariant in the leading order in temperature 0J4]]. 



Therefore, the following picture may be expected: within a correct calculational proce- 
dure the pole in the integrand of S would be somehow regularized, and the regularization 
would affect the low frequency region only. At low and medium energies the dispersion 
equation would have only one solution for the electron and only one for the positron branch. 
When the temperature grows an additional solution will appear, the same as in the high 
temperature limit. 



When calculating the Feynman diagram in quantum field theory or in the real-time 
approach to the quantum statistics a correct definition (regularization) of propagator poles 
is necessary. This is achieved by adding infinitesimal imaginary term i5 to po or p 2 depending 
on the propagator type, that makes the propagator free of singularities when 5^0. Thus a 
natural method of regularization is the use of propagators with nonzero imaginary part (this 
approach to regularization of the integrand in E was proposed in ||). The above propagator 
arises in the following way. The exact expression for the fermion self-energy is 



where S, D^, T u are the exact fermion and photon propagators and the vertex 
(within the temperature technique the 4-dimensional integral / d 4 k must be replaced by 
—^iiTJ^nl d 3 k). At high temperature, M <C |u;| <C T, where HTL approach is valid, 
the free vertex 7^ is used as T u , and propagators are free effective propagators, which are 
functions of effective masses and fermion damping (see below). At low temperature the 
self-energy must be calculated using the standard perturbation expansion. However, we 
shall replace the free fermion propagator by free damped fermion propagator, i.e. that con- 
taining OPE damping 7. For the retarded propagator it is equivalent to the substitution 
Po — ► Po + £7, thus in the temperature technique we have the standard equation with u n 
replaced by uo n + 7. For 7 at low temperature we shall take the high temperature expression 



III. ACCOUNT TAKEN OF DAMPING 




(20) 



1 



S 



7 = "^ Ins • (21) 

As a result, for £(a>) at low temperatures one obtains equations (6) - (9) again, where 
the substitution u — > ui + 27 is made (we have neglected the imaginary term 27 in the 
thermal distribution functions, since the principal contribution to the integral comes from 
the momenta k ~ T, and 7 -C A;). 

In Fig. 2 plots of imaginary and real parts of the electron branch of the self-energy 
calculated in the Feynman gauge at T = 5M, 8M and 10M are presented. One can notice 
that there are no additional solutions at positive frequency. At negative frequencies the real 
part of the inverse propagator has two zeroes, which would imply existence of two additional 
solutions. However, the right (closest to to = 0) solution to the equation 3fte(u; — M — S) = is 
always within the range of large (negative) imaginary part of £(a>), tfleE <C QmE, and thus 
the physical excitations are absent. The left zero at T = 5M and T = 8M (Fig. 2 a, b) is 
also strongly damped. However, when the temperature is increased the corresponding point 
moves leftward and at 10M escapes the high damping domain (Fig. 2 c). Therefore, 
the temperature for which an additional (propagating) mode appears is approximately 10M . 

At T ^> \u\ ^> M the fermion self-energy £ should be calculated using the effective 
electron propagator S e ff 

Q ( \ 7% (2n + l)7r 2 g 2 T 2 

Seff(p) = ~ 2 2 _ M2 , R) = * = M eff = — (22) 

Po — Pi IVI eff 1 4 

and the photon propagator D e JJ{k) 

D eff( k ) - 9»u kn ~l— U 2 ~ 9 — m\ 

U liv \ K ) ~ k 2 _ ^2^ ' K — 1 rp , t*eff ~ 12 ■ y l6 > 

In the above temperature range the self-energy in the leading order in T is 

M 2 ,, 

E M = -£L (24) 



2 In Ref. ^ the photon propagator was taken in the Aq = gauge. We are using the Feynman 
gauge, since the results are leading order gauge-invariant. 
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and the dispersion equation looks like 

M lff / ^ 

" = ^f ' (25) 
This provides the following value of OPE mass u + (energy at zero momentum) 

and CE mass oj- 

"—7=2 M <"-&2 ■ (27) 

To evaluate the localization of the additional low-frequency excitation localization one 
may use Eq. (17) with the substitution M — > M e ff ~ gT, which gives 



\UJ\ 



<ig 2 , \uj\^g 3 T (28) 



Meff 

The plot of SmE(u) presented in Fig. 3 demonstrates that the imaginary part of S is 
large in this frequency domain (the fact that insertion of OPE damping into the effective 
propagator makes 3mS(0) ~ T was mentioned in Ref. |5|]). 

Let us check whether the results for CE energy at T = 10M where the effective HTL 
scheme (a descent from high temperatures to T = 10M) is consistent with that obtained 
within the "improved" naive approach (an ascent from low temperatures to T = 10M). 
In Fig. 4 the plot of 3?eS(u;)/M e ^ (the quantity independent of temperature) calculated 
within the HTL method, and the straight line {u — M)/M e ff for T = 10 are presented. (One 
cannot yet neglect the vacuum mass in the dispersion equation at this temperature since 
M e ff and M are of same order.) It follows from Fig. 2c and Fig. 4 that the values of CE 
(and OPE) energies are perfectly matched. 

Let us emphasize that although the exact propagator has two poles one should take into 
account only one (OPE pole) when calculating S, since the residue in CE pole at k ~ T is 
exponentially small p[. 
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IV. PHYSICAL INTERPRETATION 



A satisfactory mechanism giving rise to an additional collective excitation is not revealed 
yet. In Ref. || the appearance of CE was supposed to be a result of interaction between 
the electron and a sea of electron-positron pairs. We shall show below that in the model 
describing "nonrelativistic" electron and photon, when antiparticles are absent, a CE similar 
to that in quantum electrodynamics still arises. 

Consider the Lagrangian 



L = ift{id t - e(id))il> + -(d t fd t ip - ipn 2 (id)ip) + iftij)f(iB)tp 
The self-energy S(p) of the field ip is defined as 



(29) 



S~'(p) = S^(p) + E(p) 



(30) 



where S and So are the exact and the free propagators of the field ip. At the one-loop level 
the self-energy is 



Ei(p) 



d 3 k 



f 2 (k) 



(2tt) 3 K 2 (k) - (po - e(p - k)) 2 

(the vacuum contribution is omitted). 
First, suppose that 



e(p — k) — p 1 



K(k) 



+ 



1 



(31) 



f(k) = g, e{k) = Vk 2 + M 2 , K (k) = 



(32) 



In this case Eq . (|29|) describes a model with relativistic kinematics like electrodynamics, but 
no antifermions. Let the external momentum p be zero and consider the frequency range 
M,fi <ti \u\ <C T, uj = pq. Taking into account the fact that the principal contribution 
comes from the integration over k ~ T one obtains in the leading order in T: 

g 2 T 2 



1 g 2 T 2 f°° 
Si (a?) = — —r- / xdx 
uj 4tt z Jo 



1 1 

- 1 + e x + 1 



IQuj 



(33) 



The corresponding dispersion equation may be written as 
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/ I 2 n 2 T 2 

U) = — , UJn = , 34 

uj ' 16 ' K J 

which coincides with the similar equation in QED, Eg . (p5|) . 

For negative frequencies the self-energy is a negative quantity (at M > /i), in particular, 

'e(jfe) 1 1 



2vr 2 (/i 2 - M 2 ) Jo 



(35) 



«(fc) e^ K ( fc ) - 1 e^ £ ( fc ) + 1 
(the approximate equality is valid for T ^> M, /x). Since Si (a;) is positive at high positive 
frequencies (see Eq. fl3"3|) ) and |Ei(0)| grows and becomes greater than M when the temper- 
ature increases, the situation is analogous to QED: at low temperature Ei(u;) is small and 
only one solution (corresponding to OPE) exists. When the temperature is growing two 
additional solutions appear, one of them having negative residue in the propagator pole. 
When the temperature raised further on, one of the additional solutions is shifting to large 
negative frequencies (and it would leave the damping domain if we take the OPE damping 
into account when calculating E). The frequency of the second solution remains small and 
has negative residue (the imaginary part of E would be large there if we took the damping 
into account). 

Second, consider the elect ron-phonon interaction model 

f(k) = g 2 K 2 (k)9(k D -k) , 9 2 = ^ , < k ) = ^~» > <k) = ck , (36) 

where M is the electron mass, c is the speed of sound, ko ~ Pf-, A* ~ Pf/2M, pp is Fermi 
momentum. Let the external momentum be p = pp and let us examine the frequency range 
M ^ Pp/M. In this case in the leading order in T one has 

SlM = ?£I (37, 

and the corresponding dispersion equation is 

Therefore, starting with T^pp/M the exact propagator of the field ip has poles in 
symmetric points: 
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u = ±wo = ±^T . (39) 

A typical Fermi energy is few eV, thus a pole in the symmetric point uj = —u would 
appear at 10 4 - lO 5 -^. At these temperatures the model considered is apparently 
senseless, moreover even "solids" do not exist. But we come to an important conclusion: 
just in any model, not only relativistic ones, at sufficiently high temperatures a pole in a 
symmetric point u = —u arises. That symmetric pole is a result of the growth of E and 
looks like an increase of the effective coupling constant. The latter may be explained as 
follows: when the temperature is raised the characteristic momentum grows, k ~ T, and the 
number of particles participating in the interaction (phase volume) increases drastically. 
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FIGURES 
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FIG. 2. Imaginary and real parts of £ at T = 5M, 8M and 10M (frequency and self-energy are 
normalized by fermion vacuum mass). 
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FIG. 3. 9mS as function of frequency (a; and £ are normalized by M e ff = gT/2). Within low 
frequency range \uj\/M e ff < g 2 ~ 0.01, 9mS/M e jy is large, \^smH\/M e ff ~ 5. 
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ReZ/M e ff 




FIG. 4. SKeX as function of frequency (u and £ are normalized by M e ff). Left solution 
\u>-\/M e ff 0.43 or \u>-\ = 0.64M. The latter should be compared with corresponding solu- 
tion on Fig. 2c \uj\- « 0.66M. 
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